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DARK MATTER SUBHALOS AND THE X-RAY MORPHOLOGY OF THE COMA CLUSTER 



o 

Oh- 
< 



o 

u 

43 

Ok 

6 

(N 
> 

On 

O 

m 



is 



Felipe Andrade-Santos 1,2 , Paul E. J. Nulsen 1 , Ralph P. Kraft 1 , William R. Forman 1 , 
Christine Jones 1 , Eugene Churazov 3 , and Alexey Vikhlinin 1 

1 Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138, USA 
2 Departamcnto dc Astronomia - IAG, Universidade dc Sao Paulo, Sao Paulo - SP, Brazil and 
'Max Planck Institute for Astrophysics, Karl-Schwarzschild str. 1, Garching D-85741, Germany 

Draft version April 16, 2013 

ABSTRACT 

Structure formation models predict that clusters of galaxies contain numerous massive subhalos. 
The gravity of a subhalo in a cluster compresses the surrounding intracluster gas and enhances its 
X-ray emission. We present a simple model, which treats subhalos as slow moving and gasless, for 
computing this effect. Recent weak lensing measurements by Okabe et al. have determined masses of 
~ 10 13 Mq for three mass concentrations projected within 300 kpc of the center of the Coma Cluster, 
two of which are centered on the giant elliptical galaxies NGC 4889 and NGC 4874. Adopting a 
smooth spheroidal /3-model for the gas distribution in the unperturbed cluster, we model the effect 
of these subhalos on the X-ray morphology of the Coma Cluster, comparing our results to Chandra 
and XMM-Newton X-ray data. The agreement between the models and the X-ray morphology of 
the central Coma Cluster is striking. With subhalo parameters from the lensing measurements, 
the distances of the three subhalos from the Coma Cluster midplane along our line of sight are all 
tightly constrained. Using the model to fit the subhalo masses for NGC 4889 and NGC 4874 gives 
9.1 x 10 12 Mq and 7.6 x 10 12 Mq, respectively, in good agreement with the lensing masses. These 
results lend strong support to the argument that NGC 4889 and NGC 4874 are each associated with a 
subhalo that resides near the center of the Coma Cluster. In addition to constraining the masses and 
3-d location of subhalos, the X-ray data show promise as a means of probing the structure of central 
subhalos. 

Subject headings: dark matter - galaxies: clusters: individual (A1656) - galaxies: clusters: intracluster 
medium - large-scale structure of the universe - X-rays: galaxies: clusters 



1. INTRODUCTION 

In the standard ACDM cosmology, the massive halos 
dominated by dark matter that host galaxies, groups, 
and clusters of galaxies assemble by mergers of smaller 
structures. Under the influence of gravity, uncollapsed 
matter and smaller collapsed halos fall into larger halos 
and, occasionally, halos of comparable mass merge with 
one another. Earlier generations of dark matter halos 
can survive as subhalos within the resulting aggregates, 
causing significant departures from the idealized smooth 
mass distributions often used to model dark matter halos. 
Hierarchical collapse models make clear predictions for 
the level of substructure within dark matter halos, which 
can provide tests of the ACDM model. A possible conflict 
with these predictions for galaxy sca le dark halos is the 
basis of the "missing dwarf problem" (jKlvpin et alJll99& 
iMoore et aT1ll999l : iRocha et al.ll2013| ). 

Numerical simulations of structure formation find that 
many (~ 10) dark matter halos with masses of ~ 
10 12 Mq would be accreted by a rich galaxy cluster dur- 
ing the past few Gyr. Using t h e Millennium I fe II 
simul a tions (iSpringel et ail 2005: B ovlan-Kolchin et al.l 
[2^09l) . [Fakhouri et all (120101) find that the Coma Cluster 
(Mtot = 2 x 10 15 M Q ) would have accreted ~ 7 dark 
matter halos with masses in the range 10 12 - 10 13 Mq 
in the last Gyr. Suppo rting this prediction, the analysis 
of iGiocoli et al.l (|2010f ). which identifies substructures at 
all levels of the hierarchy (subhalos. sub-subhalos, etc.), 
finds that there should be ~ 7 dark matter substructures, 
with masses in this range, present in a cluster like Coma 



today. 

Recently, iChurazov et al.l (|2012f ) analyzed X-ray sur- 
face brightness fluctuations in the central region of the 
Coma Cluster. They discuss the physical origin of these 
fluctuations, noting in particular that fluctuations in the 
gravitational potential associated with the two massive 
central galaxies in Coma could account for some of the 
larger sc ale surface brightness fluctuations in that region 
(see also I Vikhlinin et "allll994l) . In this paper, we present 
a more detailed model for the effect of subhalos, quanti- 
fied by gravitational lensing measurements, on the X-ray 
surface brightness in the central region of the Coma Clus- 
ter. 

1.1. The Coma Cluster 

The Coma Cluster (z = 0.0236) has long been consid- 
ered typical of dynamically relaxed systems. However, 
a wide range of studies have revealed that it contains 
complex substructure. Opt ical signatures of substruc - 
ture come from the work of iFitchett fc Websterl (|1987f l. 
who used a maximum likelihood method to partition the 
cluster into subsystems that have yet to c ome into dy- 
namic al equilibrium with one another, and lMellier et ail 
(Il988h . who determined mass-to-light ratios for the sub- 



structure surr ounding the two domin ant galaxies in the 
cluster core. iColless fc Dunnl (|1996| ) used a sample of 
552 galaxy redshifts to further clarify the merger his- 
tory of the Coma Cluster. X-ray studies a lso reveal 
a remarkable complexity. iBriel et ail (|1992l ) detected 
diffuse X-ray emission from the regions of NGC 4839 
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and 4911 subgroups at 6% and 1% of the total cluster 
emission, respectively. I White et al.l (|1993l) showed that 
the Coma Cluster was formed by the merging of sev- 
eral distm£t_sjibjmits_jvliich are not yet fully destroyed 
and IVikhlinin et al.1 (|1994l ) showed that the extended re- 
gions of X-ray emission in the central region of Coma 
are associated with t he subgroups NGC 4889 and 4874. 
iBiviano et al.l (|1996l ) combined ROSAT and extensive 
optical redshift and photometric samples to identify the 
main body of the cluster, suggesting that it is rotating. 
They also concluded, based on the differences in velocity 
between the NGC 4889 and NGC 4874 groups and the 
cluster mean, that these groups only arri ved recently in 
the cl uster core. In a more recent work, iGerhard et al.l 
(|2007f ) used planetary nebulae as tracers to investigate 
the ongoing subcluster merger in the Coma core. From 
the planetary nebulae line of sight velocity distribution 
they concluded that the NGC 4889 subcluster is likely to 
have fallen into Coma from the eastern A2199 filament, 
in a direction nearly in the plane of the sky, thus col- 
liding withjthe_NGC4874 subcluster coming from the 
west. lAdami et al.l (|2005fl presented a multi- wavelength 
analysis of subclustering in Coma, using X-ray data and 
a compilation of nearly 900 galaxy redshifts to partition 
the cluster into 17 subgroups that they identify as the 
remnants of infalling groups from the surrounding large- 
scale structure. 

Further evidence of ongoing infall in Coma can be 
found in X-ray observations. For exa mple, the linear 
filame nt to the southeast identified by IVikhlinin et al.l 
(|1997| ) in the Coma Cluster, which extends ~ 1 Mpc from 
the cluster center toward NGC 4911 and NGC 4921, may 
be enhanced X-ray emission due to gas stripped from an 
infalling group or the potential perturbations caused by 
tidally stripped dark matter. In either case, it is tran- 
sient on the dynamical ti me scale and so mu s t be d ue to 
recent infall. In the radio. iBrown fc Rudnickl (|2011[ ) have 
argued that a ~ 2 Mpc radio relic on the western side 
of the cluster is one signature of an infall shock, caused 
by a new, ongoing merger. If so, this merger is distinct 
from the one responsible for the formation of the X-ray 
filament. 

X-ray studies of galaxy clusters provide clues to the 
dynamical effects of local matter concentrations on the 
intracluster (ICM). The hot intracluster plasma, which 
contains most of the baryonic mass in a rich cluster like 
Coma, is perturbed by infalling groups. Regions of en- 
hanced X-ray emission define dense gas concentrations, 
for which there are a number of possible causes, includ- 
ing the gravitational effect of dark matter subhalos on 
the hot, diffuse gas. As mentioned, numerical simulations 
suggest that in the last Gyr of its history, the Coma Clus- 
ter has accreted ~ 7 dark matter halos with masses in 
the range ~ 10 12 - 10 13 M Q , which might be detectable 
via their gravitational effect on t he cluster gas. Such a n 
effect has already been detected (jVikhlinin et alJll997l ). 
However, a large-scale perturbation to the X-ray mor- 
phology has not previously been associated directly with 
embedded dark matter subhalos. Studies of the effect 
of subhalos on the gas may help us to better understand 
the cluster X-ray morphology and test predictions for the 
d istribution of dark m atter subhalos in clusters. 

lOkabe et al.l (I2010D have measured the projected mass 
distribution within 60' of the center of the Coma Clus- 



ter (around NGC 4874), using weak lensing data from 
two Subaru/Suprime-Cam fields. They detected eight 
subclump candidates, with a mean projected mass for 
severQ of the subclumps within the cluster equal to 
7.2 ± 1.9 x 10 12 M Q . In principle, analysis of the X-ray 
surface brightness of Coma can constrain the location of 
a dark matter subclump along our line of sight within 
the cluster. The gravity of a subhalo compresses the 
intracluster gas and produces an enhancement in X-ray 
brightness that depends on the density of the intraclus- 
ter gas that surrounds the subhalo. Subhalos closer to 
the cluster center will cause greater brightness enhance- 
ments simply because they are embedded in denser gas. 
To explore this possibility, we have modeled enhance- 
ments in the X-ray surface brightness resulting from adi- 
abatic compression of the intracluster gas caused by dark 
matter subhalos residing in an otherwise relaxed cluster 
potential. 

In the following sections we present the observations of 
Coma used here (Section[2]) and our model for estimating 
the enhancement in X-ray surface brightness caused by a 
dark matter subhalo embedded in the ICM (Section [3]). 
In Section |4] we present the results of numerical calcula - 
tions for three subhalos identified bv lOkabe et al.l (|2010f ) . 
followed by a discussion in Section [5] and summary in 
Section [5] The cosmology assumed here has = 0.3, 
Oa = 0.7 and Ho = 70 km s _1 Mpc _1 , implying a linear 
scale of 0.48 kpc arcsec -1 at the Coma distance of 100 
Mpc. 

2. X-RAY OBSERVATIONS AND DATA REDUCTION 

Chandra and XMM-Newton observations of the Coma 
Cluster were used make maps of the hardness-ratio and 
X-ray surface brightness for comparison with the models 
that are discussed in detail in later sections of this paper. 
In this section we outline the observations analysis. 

2.1. Chandra Observations 

The core of the Coma Cluster (A1656) has been ob- 
served several times by Chandra, for a total exposure 
time on the ACIS-I and S detectors of 120 ks (ObsIDs 
555, 556, 1086, 1112, 1113, 1114, 9714, 10672). The 
data were reduced following standard CIAO 4.4 proce- 
dures, with CALDB 4.5.3, including corrections for time 
dependence of the charge transfer inefficiency and gain. 
Data were also filtered for periods of high background 
and standard blank sky background files were used to ob- 
tain background imagcsQ Events from individual obser- 
vations were reprojected and co-added to generate com- 
bined Chandra images with a total exposure of 120 ks. 

2.2. XMM-Newton Observations 

From XMM-Newton, we used the sa me 0.5-2.5 keV 
EPIC/MOS imaging data as used by IChurazov et al.l 
(|2012t ). In brief, XMM-Newton pointings covering a field 
of more than 1° x 1° were combined into a single image, 

They quantified the contribution of background large-scale 
structure (LSS) to the projected mass distribution using Sloan 
Digital Sky Survey multi-band and photometric data, assuming 
a mass-to-light ratio appropriate for field galaxies. They found 
that one of the eight subclump candidates, which is not associated 
with any member galaxies, is significantly affected by LSS lensing. 

2 We note that the background level for the X-ray bright Coma 
Cluster is practically negligible. 
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the central part of which was fitted to our models. Peri- 
ods of background flaring were filtered out based on the 
event rate above 10 keV. Background images were made 
by scaling "blank sky fields" to match observed count 
rates in the 11-12 keV band. The image area extends 
well outside the field covered by Chandra, to the out- 
skirts of the Coma Cluster. The cleaned and calibrated 
MOS data were combined to make a raw 0.5-2.5 keV 
image, together with a background image and exposure 
map. An error map was computed assuming Poisson 
statistics for the image and background data. Simulated 
images are multiplied by the exposure map for fitting to 
the background subtracted raw image. 

3. SUBHALOS IN COMA 

As discussed in Section 11.11 simulations of hierarchi- 
cal structure formation suggest that about seven dark 
matter subhalos with masses of > 5 x 10 12 M Q are ex- 
pected to be present in the Coma Cluster now. Here we 
explore the observable impact of dark matter subhalos 
on gas near the cluster center, where their effect on the 
cluster X-ray emis sion is greatest. Using lensing data, 
IQkabe et al.l (|2010D have reported the discovery of seven 
subclumps with masses of ~ 10 13 M©, three of which 
arc marked in Figure [1] Two of these lie close to the 
giant elliptical galaxies NGC 4874 and NGC 4889, and 
their masses are consistent with the masses attributed 
to subgroups associated with these two galaxies in previ- 
ous studies (subclumps 2 and 1, respectiv ely, in Figure QJ 
iMellier et al]|1988t IVikhlinin et~alfll994[ ). We begin by 
outlining a simple model for estimating the effect of sub- 
halos on the X-ray surface brightness of a cluster, which 
we then apply to subhalos in the Coma Cluster. 

3.1. Impact on the Intracluster Gas 

We treat the subhalos as static perturbations on the 
large scale gravitational potential of the cluster. An 
off-center subhalo cannot remain at rest, but, unless 
its speed is comparable to the sound speed or greater, 
the effect of its motion on the gas is modest (e.g., 
see IMachacek et alj 120051) . The differences between the 
mean radial velocity of the Coma Cluster and the radial 
velocities of NGC 4874 and NGC 4889 are Av r = 299 
and —430 km s _1 , respectively, both small compared to 
the gas sound speed of ~ 1400 km s _1 . Although the 
transverse velocities of these galaxies are unknown, the 
lack of evidence for shock fronts or substantial asymme- 
tries in the distribution of gas in their vicinities suggests 
that their motion through the cluster is significantly sub- 
sonic. Thus, the hydrostatic approximation is likely to 
be valid for the associated subhalos. More generally, the 
motion of typical subhalos is expected to be roughly tran- 
sonic. Although transonic motion certainly alters details 
of perturbations to the gas density, the hydrostatic ap- 
proximation will provide a good estimate for the magni- 
tude of the change in the X-ray surface brightness, unless 
the speed of a subhalo is appreciably supersonic (see the 
Appendix) . 

The effect of a static subhalo is to deepen the grav- 
itational potential locally, drawing in and compressing 
the gas, thereby increasing the local X-ray emission. At 
the relevant temperatures and for typical element abun- 
dances, the broad band responses of both Chandra and 
XMM-Newton to optically thin thermal emission from 



hot gas are relatively insensitive to the gas temperature. 
For example, for a fixed emission measure, the 0.5-2.5 
keV count rate of the Chandra ACIS-I declines by a 
total of ~ 17% as the gas temperature increases from 
U-qT = 4 to 12 keV. The temperature dependence of the 
XMM-Newton response is similar. Thus, we can reason- 
ably ignore the temperature dependence of the Chandra 
and XMM-Newton responses and assume that the X-ray 
count rate per unit volume of the gas depends only on 
its density, being directly proportional to the square of 
the gas density. 

Apart fro m tiny cores of cool g as in NGC 4889 and 
NGC 4874 (|Vikhlinin et al.1 120011 ). the cooling time of 
the gas near the center of the Coma Cluster exceeds the 
Hubble time, so that the effects of radiative cooling are 
negligible. For a slowly moving subhalo, the response of 
the gas is therefore adiabatic. We assume that the per- 
turbed region is small compared to the scale over which 
gas properties vary appreciably, so that the gas affected 
by a subhalo is approximately isentropic. The momen- 
tum equation for isentropic gas may be written 



(i) 



where $ is the gravitational potential, v is the gas ve- 
locity, H = 5/cb? 1 / (2/wih) is the specific enthalpy of the 
gas, T is its temperature, /xmn is the mean mass per gas 
particle, /cb is the Boltzmann constant, and the ratio of 
specific heats is assumed to be 7 = 5/3 throughout this 
paper. For slow motion, we can ignore dv/dt to obtain 



H + $ ~ constant 



(2) 



at any given time. As discussed in the Appendix, Equa- 
tion @ provides a good estimate for the local perturba- 
tion due to a subhalo that moves at no more than tran- 
sonic speeds and has a small mass (GM t /H <C R c , the 
cluster core radius, in the notation below). The subhalo 
also needs to have been well stripped of its gas, since we 
ignore any possible effects due to the motion through the 
ICM of a remaining atmosphere bound to the subhalo. 
This is generally a good approximation for subhalos near 
the center of the Coma Cluster. Note that our assump- 
tion that the perturbations are adiabatic (isentropic) dif- 
fers from the a s sumpt ion of isothcrmality employed by 
IChurazov et al.l (|2012[ ). Isentropic gas is less compress- 
ible than isothermal gas, so that a dark halo will cause 
a smaller change in surface brightness under the model 
used here (for the same perturbation, expanded to 
first-order, the chang e in surface brightne ss estimated 
here is 60% of that of IChurazov et al.ll2012f ). 

Expressing the gravitational potential at r as ^(r) ~ 
$o( r ) + where $0 is the potential of the unper- 

turbed cluster and 5<I> is the potential of a subhalo, the 
local effect of the potential perturbation is to modify the 
specific enthalpy to H ~ Hq — (5$, where Hq is the spe- 
cific enthalpy of the unperturbed gas at the location of 
interest. For locally isentropic gas with 7 — 5/3, the 
perturbed electron density is 



n = n,Q 1 — 



(5$ 



3/2 



(3) 



where uq is the unperturbed electron density at the po- 
sition of interest. 
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a subhalo. In the absence of a magnetic field, the thermal 
conductivity would be Kg ~ 4.2 x 10 13 erg cm" 1 s" 1 K _1 . 
The effective conductivity of the magnetized ICM is 
poorly known, but expressi ng it as kq = fns , the larges t 
reasonable value is / ~ 0.2 (|Naravan &: Medvedev|[200l[ ) . 
Scaling to these values gives 



50/ 



no 



Fig. 1. — Left: positions of the three central subhalos. Subhalo 1 is centered on NGC 4889 (denoted by the number 1), subhalo 2 is 
centered on NGC 4874 (2) , and subhalo 3, which has no optical counterpart, is centered on X-ray bright source closest to the lcnsing 
position from lOkabe et all ||2010T ). Right: the XMM-Newton image contours are shown overlaid on the optical image (SAO-DSS) in the 
right panel. 

3.2. Conduction and the Adiabatic Approximation 

High thermal conductivity in the Coma ICM has the 
potential to make the perturbations more nearly isother- 
mal than adiabatic. When the effect of a subhalo can 
be regarded as a first order perturbation and the gas is 
effectively hydrostatic, first order changes in the entropy 
distribution produce no first order change in the pres- 
sure. This makes the perturbations locally isobaric in 
the presence of thermal conduction. Including conduc- 
tion in the energy equation and expanding to first order 
gives 

5p dST 2 

__ =K0 V5T, (4) 

where po, To, and kq are the unperturbed pressure, tem- 
perature, and conductivity, respectively, and ST is the 
temperature perturbation. From this, we can estimate 
the time required for thermal conduction to turn sub- 
halo perturbations from adiabatic to isothermal on the 
length scale i as 



( kT \ 5/2 / 
V8keVy V 



.1000 km s- 



0.004 cm" 3 
-l 

kpc. 



(7) 



tc 



5 Po l 2 
2kqTq 



(5) 



As a subhalo moves at speed vq through the ICM, 
the gas immediately around it is continually changing. 
The time the subhalo takes to cross a region of size I is 
Across = i/vQ. For conduction to make regions of this size 
isothermal, we must have t CO nd < Across- Equating these 
timescales determines the approximate extent, £ max , of 
the region that is kept isothermal as a subhalo moves 
through the ICM, 



2kqTq 



(6) 



Near the center of the Coma Cluster, the electron den- 
sity is close to 0.004 cm -3 (Lyskova et al., in prepara- 
tion), so that, for a temperature of 8 keV, the electron 
mean free path is ~ 3.7 kpc, which is small compared to 



For / = 0.2, this is comfortably smaller than the size 
of all three subhalos considered here at their best-fitting 
distances from the cluster center (Table [2]) . On scales 
greater than £ max , the adiabatic approximation is more 
accurate than the isothermal approximation. Because 
the long range gravitational impact of a subhalo is most 
significant for its effects on the surface brightness, the 
adiabatic approximation is appropriate for the subha- 
los considered here. Ignoring the effects of conduction 
slightly underestimates the peak surface brightness over 
a subhalo. We would need to know both the effective 
conductivity and the speed of a subhalo to correct for 
this. 

Note that the emission measure weighted temperature 
profiles in Figure [6] are quite sharply peaked and would 
be more significantly affected by a high effective thermal 
conductivity than the surface brightness profiles. Poten- 
tially, the temperature profiles can provide a probe of 
the effective conductivity on scales of a few tens of kpc, 
although existing instruments d o not have the sensitivit y 
and accuracy to use it (cf., e.g.. iMarkevitch et alj|2003() . 

3.3. Unperturbed Cluster Model 

The gas density distribution in the unperturbed clus- 
ter is assumed to have a prolate spheroidal form, with 
its major axis in the plane of the sky. With the z axis 
along our line of sight, the ellipsoidal coordinate, i?, is 
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expressed in terms of cartesian coordinates {x, y, z), with 
their origin at the cluster center, by 



x cos a 
xsmO 



R 2 =x 



+ ysmt'; 
— ycosO; 

(y' 2 + z 2 )/e 2 



where 8 determines the orientation of the major axis and 
e is related to the eccentricity, e, by e 2 = 1 — e 2 . In terms 
of R, the unperturbed gas density distribution is given 
by the /3-modcl for the electron density, 



n (R) = n c 



1+ ^ 



R 

Rc 



1 -3/3/2 



(8) 



where the central density, n c , the core radius, i? c , and j3 
are all constants. 

The temperature plays a secondary role in our mod- 
els, only entering through the unperturbed enthalpy in 
Equation ([3]). Thus, the temperature in the unperturbed 
cluster is assumed to be constant, an adequate approx- 
imation for the core region of the Coma Cluster (see 
lArnaud et aHl200H iBriel et al.ll2001[) . 

3.4. Subhalo Model 

The density distribution of gravitating matter in a sub- 
halo is assumed to have the NFW form ([Navarro et al.l 

[mi, 

= / P7T / \2 > ( 9 ) 
r / a(l + r/a) z 

where r is the distance from the center of the subhalo, 
a is the scale radius, and po is the scale density. The 
subhalo density profile is truncated at r = r t , giving a 
total mass of 



M t = 4-irpoa 3 



ln(l + w t ) - 



l + w t 



(10) 



where Wt — r t /a. The corresponding gravitational po- 
tential perturbation is 



5<Z>(r) = 



4-KGpoa 
GM t 



1 



1 + w t 



- - ln(l 

r 



r/a) 



r <r t 



r > r t . 



(11) 

The dens i ties of the subhalos determined by 
IQkabe et al.l (12010D from lensing data are higher than 
expected from simple hierarchical collapse models. Us- 
ing the NFW model, we can estimate the properties of 
a progenitor halo from those of its remnant subhalo. If 
the subhalo was simply truncated by falling into the clus- 
ter, its mass, M t , size, r t , and scale length, a, determine 
the scale density, po, of its progenitor through Equation 
(fTU|) . In terms of the mean density of the subhalo, this 
requires 



Pt 



3Mt 



Po 



3 

— <! ln(l + w) 
w a 1 



w 



w—rt I a 



--PoF{r t /a), 



(12) 



which determines po as WC U as defining the form of F(w). 
To estimate the virial radius of the progenitor halo, we 
adopt the usual approximation that the mean density 



within the virial radius is 200 times the critical density, 
evaluated at the redshift of the Coma cluster. That is, 
the virial radius of the progenitor, r v ; r , is determined by 
solving 



200 



3ff 2 
8^G 



PoF(r vil /a) = p t F(r vil /a)/F(r t /a), 

(13) 

where Equation (j!2[) has been used to eliminate po- Here 
only, H is the Hubble constant at the redshift of the 
Coma Cluster (not the specific enthalpy). The virial 
mass is then 47rp vir r^ ir /3, an increasing function of r v ; r . 
For example, for a subhalo with a mass of 10 13 Mq en- 
closed within r t = 50 kpc, if its scale radius is a = 50 
kpc, the concentration parameter, c = r V j r /a, of the pro- 
genitor group would have been ~ 19.4 (r vir ~ 970 kpc) 
and its total mass would have been ~ 1.07 x 10 14 Mq. 
Reflecting the high density of the subhalo, a concentra- 
tion parameter of 19.4 is higher than expected for such a 
massive progenitor (e.g.. lMandelbaum et al.f2008l ). Since 
the NFW scale lengths of the subhalos are largely un- 
constrained, this issue might be remedied by a different 
value for a. Differentiating Equation (|13j) with respect 
to a gives the dependence of r v i r on a as 



dlnF(w) 



d\ni 



dla i 



-r v i r / a 



dlnF{w) 



din a 
~dlnF(w) 



din; 



dln-i 



w—rt I o. 



and, from the definition of F, 



-r v i r / a 

(14) 
(15) 



where p{r) is the mean density of an NFW halo inside 
radius 7', so that 



~d\n.F{w)~ 


= 3 


>0) x 


dlnw 


w—r 1 a 


.p( r ) 



P(r v 



P(r v 



1 



d In r vir 
din a 



P(?Vir) _ p(rt) 
p(r vir ) p{ r t)' 



(16) 



Now p{r)/~p(r) < 1 and, for an NFW halo, p(r)/~p(r) is 
a decreasing function of r, so that Equation ([r5[) shows 
that r v i r is an increasing function of a. It follows that 
the virial mass determined from the subhalo properties is 
an increasing function of a. Expressing r v i r as r v i r = ca, 
where c is the concentration parameter, Equation (|16[) 
also shows that c is a decreasing function of a. The 
progenitor mass cannot be less than the subhalo mass 
Mt = 10 13 Mq and a concentration parameter of 19.4 
is already too high for that mass, so reducing a, which 
would increase the concentration parameter, cannot give 
a consistent progenitor mass and concentration parame- 
ter. Increasing a does not help cither. Raising a to 288 
kpc would make the progenitor mass 10 15 Mq, with a 
concentration parameter of 7.1, which is also too high 
for the halo mass. Since this progenitor mass is already 
about half the t otal mass of the Coma Cluster (e.g., 
iKubo et"ani2007t ). it is unreasonable for the progenitor 
to be any more massive. This leaves no reasonable pa- 
rameters for the progenitor consistent with expectations 
for standard ACDM structure formation models. 

So far, we have assumed that the remnant subhalos 
are simply truncated by falling into the cluster, but it 
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is far more likely that the mean density of a subhalo is 
reduced by its interactions with the cluster. Allowing for 
this would require even greater mean densities for the 
progenitors, exacerbating the issue of the high subhalo 
densit ies. Thus, the subhalos defined by lOkabe et al.1 
(|2010T l are smaller than expected from standard structure 
formation models, if they fell into the Coma Cluster in 
the recent past. The simplest resolution to this issue 
would be if the more massive central subhalos fell into 
Coma long in the past, when the mean de nsities of their 
proge nitors would have been higher (cf. iBiviano et al.l 
[1999 ). Our results are much less sensitive to the structure 
of a subhalo than they are to its total mass, so we do not 
pursue this issue any further here. 

3.5. Estimate for the Effect of a Subhalo 

With the approximations outlined above, the pertur- 
bation to the X-ray surface brightness must still be com- 
puted numerically (Section 2]). However, before consid- 
ering the numerical results, we make a crude estimate 
of the effect on the X-ray surface brightness of a single 
subhalo. Under our assumptions, the fractional change 
in X-ray surface brightness is 



SI 

T 



J(n 2 - n 2 )dz 
I n o dz 



(17) 



where the integrals are along the line of sight. To first 



order in (5$, Equation gives n 2 



n 2 (-35$/H ). 



For a localized perturbation, we get the further approx- 
imation, 



(n 2 — nfydz ~ Suq/Hq 



-S$dz 



6n n 



,GMt 



hi 



W 



(18) 



where w is the minimum distance of the line of sight from 
the center of the subhalo and this result applies to lines of 
sight with -no > 7*t . The integration has been truncated at 
z = ±2 max to avoid a logarithmic divergence (a spurious 
consequence of the crude approximations used here). For 
lines of sight that intersect the subhalo, i.e., for w < 
r t , the magnitude of the result is similar. Setting the 
logarithmic factor to unity and using the form ([8]) for the 
unperturbed ICM density, n , finally gives us an order 
of magnitude estimate for the fractional perturbation in 
the surface brightness near the center of the subhalo 



SI 



6GM t 



R c (R c 



7 2\3/3-l/2 



2)30 



I R c H B(3/3 - 1/2, 1/2) (R 2 + z 2 ) 

(19) 

where B(a, b) is the f3 function, and vo v and z p are cylin- 
drical coordinates for the center of the subhalo measured 
with respect to the cluster center, i.e., their projected 
separation and their separation along our line of sight, 
respectively (Figure [2]). For an unperturbed gas temper- 
ature of fceTo = 8 keV, a subhalo mass of M t = 10 13 M , 
and /3-model core radius of R c = 300 kpc, the first fac- 
tor on the right in Equation ([H]) is ~ 0.17 for f3 ~ 2/3. 
Thus, such a subhalo should be detectable in X-ray ob- 
servations, unless it is far from the cluster center (|z p | or 
tup ^> R c ), Note that the estimate for SI /I in Equation 



(119)) only depends on the mass of a subhalo, not on its 
size or scale parameters. This reflects the insensitivity of 
the full, numerical results to details of the mass distri- 
bution in the subhalos. It is a consequence of the long 
range of the gravitational force. 

4. SIMULATED X-RAY IMAGES OF THE COMA CLUSTER 

In this section we present numerical results, employing 
the model outlined above to compute the effects of subha- 
los on the X-ray surface brightness of the Coma Cluster. 
The X-ray surface brightness is computed by integrat- 
ing the square of the gas density, n, given by Equation 
((3]), along lines of sight through a model cluster with 
an unperturbed gas temperature of fee To = 8 keV and 
an unperturbed gas density profile given by the /3-model 
(Section I3.3[) . Most parameters of the /3-model are de- 
termined in the fits, as outlined below. Gas densities are 
evaluated and integrated numerically on a rectangular 
parallelepiped of 321 x 321 x 1000 cubic cells, 1.906 kpc 
on a side and centered on the cluster. The long axis of 
the grid is parallel to the line of sight. 

Subhalos are modeled with NFW p r ofiles (Section l3. 41) . 
using parameters from lOkabe et al.l ( 2010f ) . In the fol- 
lowing analysis, apart from Section 14.21 the NFW scale 
radius, a, is set to 50 kpc in a subhalo that is also trun- 
cated at r t = 50 kpc to be consistent with the lensing 
results. The magnitude of the change in surface bright- 
ness due to a subhalo increases with the density of the 
surrounding gas, so that it is maximized when the sub- 
halo is closest to the cluster midplane, defined by z p = 0. 
Because the effect of a subhalo is sensitive to its distance 
from the midplane, we can constrain the locations of the 
subhalos along our line of sight. 

4.1. Fitting 

Models are multiplied by the exposure map and added 
to the background image for fitting to the raw 0.5- 
2.5 keV XMM-Newton image of the central part of the 
Coma Cluster, made from the combined EPIC MOS 
data. Goodness of fit is determined using x 2 as the fit 
statistic. The image was binned into 2" x 2" pixels and 
the region fitted had an extent of 10'.7 x 10'.7. The error 
estimate for each pixel includes Poisson errors for the raw 
image and for the background data. The fit was masked 
to exclude point sources evident in the XMM-Newton 
image (Figure [2]). Away from point sources, the average 
photon count per pixel in the raw image varies from over 
50 near the center to about 8 at the periphery. 

Models include the three subhalos marked in Figure [TJ 
The centers of the subhalos corresponding to the giant 
elliptical (gE) galaxies NGC 4889 and NGC 4874 were 
placed on the sky at their respective optical centers, 
rather than at the locat i ons de termined from the lens- 
ing data of lOkabe et al.1 (|2010T l . The center of subhalo 
3 was placed at the location of the X-ray bright source 
closest to the lensing position, since the X-ray surface 
brightness is elongated toward this position. The lensing 
positions are consistent with these positions, within the 
errors of the lensing measurements, but using the latter 
positions gives somewhat better fits (A% 2 ~ 700, Ta- 
ble^. Coordinates used for the subhalo centers are, for 
NGC 4889 (J2000) = (195.034°, 27.977°), for NGC 4874, 
(194.899°, 27.959°), and for the southern mass concen- 
tration, (195.006°, 27.855°). Fixed parameters for the 
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Cluster 
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Fig. 2. — Left: parameters defining the location of a subhalo. Right: XMM-Newton image showing the mask used to remove the point 
sources. 



TABLE 1 

Parameters of Subhalo Models 



Subhalo 


a J2000 


<5j2000 


Aft 


a 


n 




(deg) 


(deg) 


(10 12 M ) 


(kpc) 


(kpc) 


1, NGC 4889 


195.034 


27.977 


11.0 


50 


50 


2, NGC 4874 


194.899 


27.959 


6.57 


50 


50 


3 


195.006 


27.855 


7.57 


50 


50 



subhalo models are given in Table [T] The location of 
the center of the X-ray emission from the Coma Clus- 
ter was determined by fitting a circularly symmetric f3- 
modcl to the central region of the cluster. This center, 
(194.932°, 27.929°), was kept fixed in all subsequent fits. 
All other parameters for the cluster /3-modcl (9, e, R c , 
P, normalization) are determined in the fits. For each 
subhalo, the only free parameter is its location along our 
line of sight, measured by its distance from the cluster 
midplane (z p in Figure [2]). 

Fitting models with up to eight nonlinear, free pa- 
rameters with ~ 10 5 degrees of freedom (dof) presented 
some challenges. To ensure that the global best fit 
was found, two differe nt fitting procedures, a Leven berg- 
Marquardt algorithm ([Bevington fc Robin son 2003]) and 
a grid search, were employed to cross check one another. 
Initially, the three subhalos were placed at z = 1 Mpc, ef- 
fectively removing them from the fit, and the best fitting 
parameters for the /3-modcl alone were determined. This 
fit gives P = 0.80, R c = 351 kpc, e = 0.83, 9 = -3.9°, 
and x 2 = 139204.6, or a reduced x 2 (x 2 /dof) of 1.3881 
(Table O. 

When the three subhalos are included in the fit, with 
subhalos 1 and 2 at the locations of the two gE galax- 
ies, allowing the positions along the line of sight of the 
three subhalos, z,; for i = 1,2,3, to vary simultane- 
ously gives the best fitting parameters z\ = 121 kpc, 
z 2 = 7.9 kpc, z 3 = 399 kpc, R c = 387 kpc, 9 = -0.1°, 
e = 0.84, for which x 2 = 119114.7, making the reduced 
chi squared x 2 /dof = 1.1878. The improvement in the 



fit, of Ax 2 = 20089.9 for the addition of three parame- 
ters, is highly significant. Visual comparison of the model 
and cluster images in Figure [3] shows that the subha- 
los, particularly the two associated with NGC 4889 and 
NGC 4874, make the model X-ray image very similar to 
the observed cluster. Using A% 2 = 2.706 to estimate 
confidence ranges for the three subhalo positions, with 
all of the other fitted parameters free, gives 90% confi- 
dence ranges of z\ = 121 ± 5 kpc, z 2 = 7.9 ± 7.5 kpc, and 
z 3 = 399 ± 10 kpc. 

4.2. Structure and Masses of the Subhalos 

As discussed in Section 13.41 the subhalo properties 
determined from the lensing measurements give them 
higher densities than expected for hierarchical collapse 
models. To test how sensitive our results are to the struc- 
ture of the subhalos, we varied the truncation radius of 
the most massive subhalo, corresponding to NGC 4889 
(subhalo 1), while its mass was kept fixed. Since the 
NFW model is unlikely to be valid for a truncation ra- 
dius, r t , smaller than the NFW scale length, a, we tied 
these two parameters together for these fits (a — r t ). All 
parameters for the other subhalos were frozen at their 
best-fitting values, since the fits are not sensitive to these. 
The resulting x 2 and positions along our line of sight for 
subhalo 1 are plotted against r t in Figure |4j The fit is 
clearly affected by the structure of this subhalo, with a 
moderately significant improvement of Ax 2 — 4 for the 
best fit. However, at 47 kpc, the best fitting truncation 
radius for subhalo 1 is actually smaller than the value 
of 50 kpc from the lensing results, only exacerbating the 
problems caused by the high density of this subhalo (Sec- 
tion 13. 4j) . This result depends on our model through 
properties that are not well determined, so it should be 
treated with caution. For example, Table [2] shows that 
the best fitting position, z\, for the NGC 4889 subhalo is 
quite sensitive to the position on the sky of this subhalo. 
Nevertheless, the fact that the subhalo structure has an 
appreciable effect on the fit shows that the X-ray data 
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TABLE 2 

Parameters of X-Ray Surface Brightness Model 



Model" 





c 


e b 


Rc 


Zl 


Z2 


23 


dof 


x 2 








(deg) 


(kpc) 


(kpc) 


(kpc) 


(kpc) 






Cluster only c 


0.80 


0.83 


-3.9 


351 


1000 


1000 


1000 


100283 


139204.6 


Subhalos (lensing) £i 


0.82 


0.84 


3.2 


395 


39 


5.8 


361 


100280 


119817.0 


Subhalos (gE) e 


0.80 


0.84 


-0.1 


387 


121 


7.9 


399 


100280 


119114.7 



Note. — "Using the cluster /3-model of Section 13.31 and subhalo models of Section 13.41 
'Position angle of the major axis, increasing north of west. c Subhalos fixed at 1 Mpc from 
the cluster midplane. d Subhalos positions on our line of sight free to vary, with the subhalos 
centered at the lensing positions. e As for d, but with subhalos 1 and 2 centered at the positions 
of the two gE's and subhalo 3 centered on the X-ray bright source. 



200 kpc 




Fig. 3. — Left: Model X-ray image for the Coma Cluster, made by embedding three dark matter subhalos in the elliptical /3-model of 
Section 13.31 Right: 0.5-4 keV XMM-Newton X-ray image of the Coma Cluster. North is up and east to the left. Increased brightness due 
to the presence of the subhalos elongates the model image in the direction of each subhalo, particularly toward the east. These features 
correspond well to structure in the X-ray image. 



could potentially be used to probe its structure. Results 
for the NGC 4874 subhalo are similar. 

From Section 13. 5[ the X-ray surface brightness of the 
model cluster is expected to be considerably more sensi- 
tive to the total mass of a subhalo than to its structure. 
In principle, we can determine the subhalo masses by 
fitting the X-ray image, using data that is completely 
independent of the lensing mass measurements. To in- 
vestigate the agreement between the lensing and X-ray 
mass determinations, we have again used subhalo models 
with a = r t allowed to vary, but here the scale density, 
Po, or, equivalent ly, the mean density of the subhalo, is 
kept fixed, so that its mass varies as oc . Parameters 
for the cluster model are free to vary, but those for the 
other subhalos are kept fixed at their best-fitting values 
(last row of Table ©. 

The left panel of Figure [5] shows results for varying 
the mass of the NGC 4889 subhalo (subhalo 1). Plot- 
ted against r t , the subhalo mass is given in the bottom 
panel, distance from the cluster midplane in the middle 



panel, and x 2 m the top panel. Reducing the mass of the 
subhalo decreases its impact on the gas, hence the X-ray 
surface brightness, but this can be offset by shifting the 
subhalo toward the Coma Cluster midplane, where the 
unperturbed gas density is greater. The best-fitting mass 
for subhalo 1 is 9.14 x 10 12 M w , c± 1.2a lower than the 
lensing mass of IQkabe etHI (|2010D . The corresponding 
improvement in the fit is A% 2 ~ 237. Results for vary- 
ing the mass of the NGC 4874 subhalo are shown in the 
right panel of Figure El Here, the best fitting mass of 
7.61 x 10 12 M Q is ~ 0.9cr greater than the lensing mass 
and the improvement in the fit is A% 2 ~ 242 (sec Table 
[3]). This subhalo starts from much closer to the cluster 
midplane than subhalo 1 and the nonmonotonic variation 
of Z2 seen in the plot results from the interplay between 
the subhalo and cluster parameters (9 and e in particu- 
lar, see Section |3~3"|) . The mass of subhalo 3 also affects 
the fit significantly, but \ 2 continues to improve as it is 
increased well beyond the range consistent with the lens- 
ing data (e.g., doubling the subhalo mass moved it to 
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Z3 ~ 510 kpc from the midplane and an improvement of 
A\ 2 — 350). This indicates that our model assumptions 
are poorer for this subhalo. 




Fig. 4. — Effect of varying the truncation radius of subhalo 1. 
This plot shows how the fit depends on the structure of subhalo 1, 
through its truncation radius, rt. The lower panel shows the best 
fitting distance of the subhalo from the cluster midplane along our 
line of sight, while the upper panel shows x 2 ~ 119000. For these 
fits, the scale radius of the subhalo, a, is tied to its truncation 
radius, rt, the subhalo mass is fixed and the parameters of the 
other subhalos are frozen. The remaining fit parameters are free. 



4.3. Temperature Perturbations 

Next to the gravitational effect of subhalos, the most 
likely cause of large-scale enhancements in the X-ray sur- 
face brightness of the Coma Cluster is low entropy gas 
that has fallen in during the ongoing collapse. In that 
case, the pressure of the gas responsible for the excess 
emission would be similar to that of its surroundings, but 
its density would be higher, boosting the X-ray bright- 
ness. By contrast, gas compressed adiabatically in the 
potential of a subhalo should be hotter than its sur- 
roundings. For the model used here, the temperature 
enhancement due to adiabatic compression can be deter- 
mined from H = H a — <5$ (Equation ^f). Using the 
emission measure weighted temperature to estimate the 
fractional temperature perturbation, the model gives 



ST 
T 



Jn 2 Hdz 
Hq J n 2 dz 



1 



fnl(-6$/H )dz 
Jn^dz 



where the third form is obtained by expanding to first 
order in 5$/ Hq, using Equation (j3|), and the last form is 
obtained in the same manner from the expression for the 
surface brightness fluctuations in Equation (fTT)) . Com- 
bining this result with the rough estimate for SI /I from 
Section 13.51 the maximum temperature fluctuation due 
to a subhalo is expected to be roughly 5%. Evaluated nu- 
merically for the best-fitting parameters in the last line of 
TableEl the maximum temperature rise is only about 3% 
for NGC 4889 and less for the other two subhalos (Figure 
[5]) . These are not large enough to detect in small regions 
with the current data, particularly at the relatively high 
temperature of the Coma Cluster. 



However, we can check for signatures of cooler gas that 
might be responsible for the enhanced X-ray emission. 
For this purpose, we have made maps of the hardness 
ratio in 50 kpc regions around each of the subhalo loca- 
tions by binning Chandra data for the bands 0.3-1.5 keV 
and 1.5-6.0 keV into 15."74 ~ 7.5 kpc pixels. The im- 
ages in Figure [7] show (hard — soft)/ (hard + soft). Apart 
from the innermost pixels for NGC 4889 and NGC 4874, 
ther e is no indication of co oler gas in any of the subha- 
los. iVikhlinin et al.1 (|2001l ) showed that there are ~ 3 
kpc cores of 1-2 keV gas in both of these gE galaxies, 
which affect the central pixel of the respective hardness- 
ratio images. These regions are masked out of the image 
used to fit the models and they certainly do not account 
for the extended regions of enhanced X-ray emission as- 
sociated with the subhalos. There are no other signs 
of cooler gas in the ima ges of Figure Cp consistent with 
previous investigations flArnaud et al.l l200lt iBriel et all 
120011 : iNeumann et al.ll2003|) . In particular, there is no 
sign of a radial temperature increase in any of these re- 
gions. Apart from the solitary low pixels for the two 
cool cores, the histograms of the hardness ratio for each 
of these regions are reasonably consistent with a con- 
stant value. This is illustrated in Figure 13 which shows 
histograms of the hardness ratio for the regions around 
NGC 4889 (left panel and panel b) and NGC 4874 (right 
panel and panel c) in Figure [7] Statistical noise in the 
soft and hard counts contributes to the standard devia- 
tion of the hardness ratio as 



CHR,stat 



2 v /7^ 



(h 



(21) 



where s and h are the mean counts for the soft and hard 
bands, respectively, and a s and c/j are the respective 
standard deviations for the two bands. Dashed lines in 
Figure [S] show Gaussians with means, ^hr, equal to the 
data mean and normalized to match the number of sam- 
ples with standard deviations given by Equation (|21[) . 
The plotted number of samples per bin is given by 



A 7 . 



Nw 



samples 



CHR,statv27r 



cxp 



(ip - Mhr) 2 

zo HR,stat 



(22) 



where ip is the hardness ratio, N is the total number of 
samples, and w is the width of the histogram bins. The 
agreement between these Gaussians and the histograms 
shows that, apart from the single cool pixel in each core, 
the histograms show little evidence for variation of the 
hardness ratio over each subhalo. 

5. DISCUSSION 

Under our model assumptions, the distances from the 
cluster midplane of each of the three central subhalos are 
well constrained at the 90% confidence level (zi = 121±5 
kpc, z 2 = 7.9 ± 7.5 kpc and z% = 399 ± 10 kpc), provid- 
ing useful information on their three-dimensional loca- 
tions within the Coma Cluster (Section 14. 1|) . In par- 
ticular, subhalo 1 (NGC 4889) is closer to the cluster 
midplane than its distance from the cluster center and 
subhalo 2 (NGC 4874) is even closer to the cluster mid- 
plane. By contrast, subhalo 3 is 1.7 times further from 
the midplane than it is from the cluster center. These 
results arc affected by our model assumptions, particu- 
larly those concerning the masses and locations on the 
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Fig. 5. — Effect on the fit of varying subhalo mass. Left: Fit results for the NGC 4889 subhalo are plotted against its truncation radius, 
rt, with the mean density of the subhalo held constant. Right: same for NGC 4874. The bottom panels show the subhalo mass (oc rj? ), 
the middle panels show its position along our line of sight, and the top panels show x 2 ~ 119000. Model parameters match the lcnsing 
parameters in Table [T1 for rt = 50 kpc. For these fits, the scale radius of the subhalo, a, is tied to its truncation radius, rt, the mean density 
of the subhalo is fixed, while the parameters of the other subhalos are frozen. 



TABLE 3 
Best-fit Masses of the Subhalos 



Subhalo «J2000 feooo Mt a rt 

(deg) (deg) (10 12 Mq) (kpc) (kpc) 



1, NGC 4889 195.034 27.977 9.1 47 47 237 

2, NGC 4874 194.899 27.959 7.6 52.5 52.5 242 



i 1 r 
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Fig. 6. — Emission measure weighted temperature profiles com- 
puted from our model for the three subhalos. The peak temper- 
ature over the NGC 4889 subhalo is 3% greater than the model 
temperature of 8 keV. 

sky of the subhalos. They are also affected by our largely 
arbitrary assumption that the gas distribution in the un- 
perturbed cluster is prolate, with its major axis in the 
plane of the sky. While these assumptions affect the dis- 
tances of the subhalos from the cluster midplane, they 
do not alter the qualitative result that subhalos 1 and 
2 lie close to the cluster midplane, while subhalo 3 does 
not. Our results are consistent with previous findings 
that NGC 4889 and NGC 487 4 lie close to the cluster 
center (e.g. iBiviano et ai1ll996l ). 

Nonlinearity, from Equations (|3|) and (Ti~7)) , means that 
the impact of two or more subhalos on the X-ray sur- 



face brightness increases when they are closer together. 
As a result, the placement of subhalos about the clus- 
ter midplane can alter the fit. This is seen in the fit 
results for subhalos 1 and 2. In computing the confi- 
dence range for subhalo 2, we found that A% 2 continues 
to rise as z-i decreases through zero. In principle, the fit 
should be symmetric about z-i = 0, but, to enforce that 
result it would be necessary to flip the positions of the 
other subhalos to the opposite side of the cluster mid- 
plane as subhalo 2 crossed the midplane. This behavior 
shows that the fit is better (formally, at the 90% level) 
when subhalos 1 and 2 lie on the same side of the cluster 
midplane. Of course, this does not determine which side 
of the midplane that is. Subhalo 3 is too far from the 
other two subhalos to obtain any information about its 
location relative to them. 

When the masses of the subhalos were allowed to vary 
(Section 14. 2[) . the best fitting masses for the subhalos 
associated with the two gE galaxies near the cluster cen- 
ter (mi = 9.1 x 10 12 M Q and m 2 = 7.6 x 10 12 M Q ) 
were found to be consistent with their lcnsing masses 
(m Uons = (11.0 ± 1.5) x 10 12 M e and m 2 , lcns = (6.6 ± 
1.2) x 10 12 M Q ), though not for subhalo 3. Formally, the 
large improvements in the fit show that the X-ray masses 
determined for the NGC 4889 and NGC 4874 subhalos 
would be much more accurate than the lensing masses. 
However, given the shortcomings of our model, particu- 
larly the approximation that the subhalos are static, the 
results would be dominated by systematic uncertainties. 
More accurate and more realistic models are required to 
assess the systematics, so we make no attempt to pursue 
this further here. The remarkable consistency between 
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Fig. 7. — Hardness-ratio maps for the subhalos. a) a representative region lying between the 3 subhalos. b) subhalo 1, centered on 
NGC 4889. c) subhalo 2, centered on NGC 4874. d) subhalo 3, centered on the X-ray bright source. The hardness ratio is defined in the 
text. Each map covers a region 102" or ~ 50 kpc in radius. Apart from the single cool central pixels for NGC 4889 and NGC 4874, there 
is no significant sign of a radial temperature gradient in any of these regions. 
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Fig. 8. — Histograms of the hardness ratio for regions b (NGC 4889, left panel) and c (NGC 4874, right panel) of Figure[7] The dashed 
lines show Gaussians with the same means as the data, but with standard deviations calculated to include only the effect of Poisson noise 
on the hardness ratio. Apart from a single cool central pixel at the left in each histogram, they are reasonably consistent with a constant 
hardness ratio in each region. 

the X-ray and lensing masses for the two central subhalos 
is noteworthy. It lends strong support to the interpreta- 
tion of these mass concentrations as subhalos associated 
with the two gE galaxies, which lie close to the midplanc 
of the Coma Cluster. 

The lack of any significant amount of cool gas as- 
sociated with the subhalos is consistent with these re- 
sults. However, more sensitive observations are needed 
to detect the temperature enhancements that should 
be caused by adiabatic compression associated with a 
subhalo. With existing X-ray instruments, the small 
temperature rise should be more readily detectable in 
cooler clusters and groups. Apart from providing further 
constraints on subhalo properties, detecting a tempera- 
ture rise will help to distinguish gasless subhalos from 
other sources of X-r ay enhancements in clusters (e.g. 
iChurazov et al.ll2012D . 



6. CONCLUSIONS 

We have presented an approximate model for calculat- 
ing the impact of initially gasless subhalos on the X-ray 
image of a cluster which relies on the subhalos being slow 
moving and the gas response being adiabatic. The model 
was applied to XMM-Newton and Chandra data for the 
Coma Cluster, using properties determined from lensing 



data fo r three central subhalos reported by lOkabe et al.l 
(|2010f ). Two of the subhalos correspond to the gE galax- 
ies, NGC 4889 and NGC 4874. When combined with an 
ellipsoidal model for the gas distribution in the unper- 
turbed cluster, including the effects of these three sub- 
halos produces a model X-ray image that is a much bet- 
ter fit to the observed cluster than the ellipsoidal cluster 
model alone. 

Using the lensing masses for the three subhalos, the 
distance along our line of sight of each subhalo from the 
cluster midplane is tightly constrained, with 90% confi- 
dence ranges no more than ±10 kpc (zi = 121 ± 5 kpc, 
z 2 = 7.9 ± 7.5 kpc, and z 3 = 399 ± 10 kpc). The fits 
are improved considerably when the subhalos associated 
with the two gE galaxies are assumed to be centered on 
those galaxies, rather than at the lensing positions, which 
are consistent within the errors of the lensing determina- 
tions. The model favors the subhalos associated with the 
two gE's lying on the same side of the cluster midplanc, 
further constraining their locations in three dimensions. 
The results also show the potential for using X-ray im- 
ages to constrain the internal structure of the two central 
subhalos. 

For each of the two gE galaxies, freeing the total mass 
of the subhalo gave a significantly better fit (A% 2 ~ 
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240). For the NGC 4889 subhalo, the best fitting mass 
is 9.1 x 10 12 M Q (~ 1. 2 standard deviati ons smaller 
than the lensing mass of lOkabe et "all 120101 ) . while, for 
NGC 4874, the best fit is 7.6 x 10 12 M (~ 0.9 standard 
deviations more than the lensing mass). Formally, the 
X-ray masses are much better determined, although the 
results are dominated by systematic uncertainties that 
can only be circumvented with more accurate models. 
The X-ray and lensing mass determinations rely on com- 
pletely independent data, so that the agreement between 
them supports the assumptions that underlie both, most 
critically that these two subhalos reside near the center 
of the Coma Cluster. Fitting the X-ray image does not 
determine a mass for the third subhalo, suggesting that 
our model assumptions may be less accurate for it. 



Our results highlight the potential for using high qual- 
ity X-ray data to probe substructure in galaxy clusters. 
Particularly when used in combination with lensing mea- 
surements to constrain the mass of a subhalo, X-ray data 
can constrain where along our line of sight a subhalo lies 
within a cluster. The X-ray data can also be used to de- 
termine subhalo masses and, potentially, to probe their 
internal structure, i.e., their truncation and scale radii. 
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APPENDIX 
LIMITS OF THE APPROXIMATIONS 

Two conditions must be satisfied in order for the approximations used here to be valid. First, gas displacements 
caused by a subhalo need to be small compared to t he sc ale of entropy variations in the unperturbed atmosphere for 
the perturbed gas to be treated as isentropic (Section l3.1[) . Second, for the fluid acceleration, dv/dt of Equation (p]), to 
be negligible, it must be small compared to the gravitational acceleration due to the subhalo, V(5<3?. Here, we determine 
conditions for these approximations to hold. 

The first condition will be satisfied if displacements are small compared to the core radius, i? c , of the unperturbed 
atmosphere (Equation ([8])). We can estimate the displacements caused by a slowly moving subhalo embedded in 
isentropic gas from the requirements of mass conservation. Assuming that displacements around a subhalo are radial 
with respect to the center of the subhalo, the mass of gas within any sphere fixed to the fluid elements of the gas is 
fixed. In the notation of Section [3. 11 conservation of mass therefore requires 

^-n r 3 =[ n{r')Anr' 2 dr\ (Al) 
Jo 
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where the radial coordinates, r and r' , are centered on the subhalo and Sr is the inward radial displacement of the 
gas. Using Equation ([3]) and expanding to first order in the perturbation, this gives 

\n [r 3 - (r - Sr) 3 } ~ n Q r 2 Sr / [-S<P{r')}r' 2 dr' . (A2) 

Outside the truncation radius, r t , the gravitational potential for our subhalo model ( Section I3.4[) is simply (5$(r) = 
—GM t /r, with M t constant. Although this form overestimates \5&\, hence Sr, for r < r t , the spherical geometry 
ensures that the difference it makes to our estimate for Sr diminishes rapidly with increasing radius outside r t . Using 
(5$(r) = —GM t /r in Equation (|A2[) to estimate Sr (an overestimate), gives 

independent of r. Although the displacements around a moving halo will not be purely radial, this provides a good 
estimate for the magnitude of the actual displacements. For the subhalos considered here, GM t /Ho <§; R c and the 
first requirement is well satisfied. 

To estimate the fluid accelerations due to a moving subhalo, consider a frame in which the subhalo remains at rest, 
centered on the origin. At large distances from center of the subhalo, the fluid moves through this frame with velocity 
Vo, where — vo is the velocity of the subhalo with respect to the cluster. In the absence of the subhalo, the path of 
a fluid element would be given as a function of th e tim e by ro(t) = rj + vot, where ri is a constant. Including the 
approximate radial displacement, Sr, of Equation (|A3[) , in the presence of a subhalo the path of the fluid element is 
given approximately by 

ri (t) = r (t)-Sr^§-. (A4) 
|ro (tj| 

Differentiating this with respect to the time once gives the fluid velocity 

dri(t) Sr , . , , , . 

~dt " v ° ~ mop x x (A5) 

and twice gives the acceleration 

d ±i (t) St St 

dt 2 = 3 |r (t)| 5 V ° ' r °^ r °W X ^ Vo X ~ | ro ( t )|3 v ° x t v ° x r °(*)]- ( A6 ) 

In terms of the unit vectors Vo = vo/|vo| and ro{t) = ro (i) /|ro (i) | , the cosine of the angle between Vo and ro(t) is 
it = vo ■ ro(t). We then have |ro(i)| 2 = b 2 /(l — u 2 ), where b is the distance of closest approach of i"o(i) to the center of 
the dark matter halo (the impact parameter of the fluid clement). With this notation, the acceleration can be written 

as 

d -^£ L = JK| 2 (1 - u 2 ) [2 U v„ + (1 - 3u 2 )r (t)] • (A7) 
Its magnitude is maximized when u = 0, i.e., when the fluid is closest to the center of the subhalo, giving 

d 2 r! Sr\v \ 2 3GM t |v | 2 



dt 2 



b 2 4H Q b 2 



(A8) 



The gravitational acceleration due to the subhalo at this point is approximately GM t /b 2 , so the requirement that 
the fluid acceleration be small compared to the gravitational acceleration is simply 

|v | 2 «^ = 2 s 2, (A9) 

where we have expressed the specific enthalpy in terms of the squared sound speed (Hq = 3sq/2 for 7 = 5/3). This 
requirement imposes the condition that the speed of the subhalo should be transonic at most. 



